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Abstract 

We study how quantum fluctuations of the metric in covariant Horava-Lifshitz grav¬ 
ity influence the propagation of classical fields (complex scalar and photon). The effective 
Lorentz-symmetry violation induced by the breaking of 4-dimensional diffeomorphism is then 
evaluated, by comparing the dressed dispersion relations for both external fields. The con¬ 
straint of vanishing 3-dimensional Ricci scalar is imposed in the path integral, which there¬ 
fore explicitly depends on two propagating gravitational degrees of freedom only. Because 
the matter fields are classical, the present model contains only logarithmic divergences. Fur¬ 
thermore, it imposes the characteristic Horava-Lifshitz scale to be smaller than 10 10 GeV, if 
one wishes not to violate the current bounds on Lorentz symmetry violation. 


1 Introduction 

Among the different approaches to Modified Gravity, Horava-Lifshitz gravity (HL) con¬ 
sists in choosing different scaling properties of time and space coordinates, and therefore 
breaking 4-dimensional diffeomorphism |lj. This feature allows to increase the number of 
space derivatives of the metric, while keeping the number of time derivatives minimal. This 
leads to a better ultraviolet (UV) behaviour of graphs, without introducing ghosts, and 
therefore provides a power-counting renormalisable theory of gravity. 

A problem of the original model of HL gravity is the existence of an additional scalar 
degree of freedom (dof) for the metric, which can be understood as a Goldstone mode arising 
from breaking of 4-dimensional diffeomorphism [2]. A solution to this problem has been 
proposed in [3], where an auxiliary held A is introduced, such that its “equation of motion” 
leads to an additional constraint, eliminating the unwanted scalar degree of freedom of the 
metric. The resulting theory is invariant under a new Abelian gauge symmetry 17(1) which 
involves the metric components, the auxiliary held A and an additional auxiliary held. This 
gauge symmetry, together with the 3-dimensional diffeomorphism of HL gravity, can be 
shown to be equivalent to a 4-dimensional diffeomorphism at the lowest order in a post- 
Newtonian expansion, showing the equivalence with GR at long distances. For this reason 
this extension is called covariant HL gravity. We note however, that the long-distance limit 
is not obviously recovered: it has been shown in [4] that the equivalence principle is not 
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automatically retrieved in the infrared, and that the meaning of the auxiliary fields and their 
coupling to matter are still open questions. 

Because of the anisotropy between space and time, HL gravity is naturally described in 
terms of the Arnowitt-Deser-Misner (ADM) decomposition of the metric, which expresses a 
space-time foliation. An important consequences of space-time anisotropy is the possibility of 
imposing the lapse function N to depend on time only. This situation is called the projectable 
case and will be considered in this article. The non-projectable case has been studied in p£j 
for static spherically symmetric solutions of covariant HL gravity. 

This covariant extension to HL has led to several studies, including spherically symmet¬ 
ric solutions |6j and their relevance to an alternative model for galaxy rotation curves [7], 
cosmological solutions [8], as well as theoretical and phenomenological consistency tests of 
the theory [9]. 

In the present article, we study the effective dispersion relation for classical external fields 
(scalar and photon), coupled to the covariant extension of HL gravity, after integrating over 
gravity degrees of freedom at one-loop, on a flat background. The constraint provided by 
the auxiliary field A consists in setting the curvature tensor to zero, R = 0, which we impose 
in the path integral over metric fluctuations. The restoring force for quantum fluctuations 
is then provided by higher order space derivatives of the metric, and leads to a prediction of 
the Lorentz-symmetry violating effective dispersion relations for these fields. 

Studies of effective dispersion relations for Lifshitz-type models in flat space time have 
shown the limitations of phenomenological viability of these models. It was first noted in 
PI. where the unnatural fine tuning of bare parameters is shown, in order to match the 
light cones seen by two different scalar particles interacting. Similar studies were done in 
mi, where the effective dispersion relation for interacting Lifshitz fermions is derived, in the 
case where flavour symmetry is broken. Similarly, [L2] describes fermion effective dispersion 
relations in Lifshitz-type extensions of Quantum Electrodynamics. 

The effective speed of light seen by matter interacting with HL gravity is studied in [ T3j . 
The Authors derive this effective speed seen by a scalar field and an Abelian gauge field, and 
compare these to measure Lorentz-symmetry violation. We are doing a similar study here, 
treating matter as classical though, and imposing the constraint R — 0 in the integration 
over graviton dof. We note here that R = 0 is not a gauge choice, but a constraint which 
has the physical effect to remove one degree of freedom in the theory. As described below, 
the constraint leads to the vanishing of one of the scalar components of the space metric. 

We note here related works, involving quantization of HL gravity. [14] uses Exact Renor¬ 
malisation Group methods to provide an insight into the existence of asymptotically safe 
gravity [15]. Also, [16] is based on causal dynamical triangulations to study phase transi¬ 
tions of space time geometry in (2+l)-dimensional HL gravity. Finally, [lTj describes how 
(l+l)-dimensional HL gravity can be quantised in a similar way as a harmonic oscillator. 

A short review on the covariant version of HL gravity is presented in the next section, 
and the coupling to matter fields is presented in section 3, together with the integration of 
gravitons. This calculation is done for A ^ 1 (parameter in the kinetic term of gravitons), but 
the result does not depend on A, which is a consequence of the vanishing of the trace of the 
fluctuating space metric. Section 4 presents a phenomenological analysis and our conclusions. 
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2 Covariant Horava-Lifshitz gravity 

2.1 The action 

We consider the anisotropic scaling z — 3, such that the mass dimensions of time and 
space are [£] = —3 and [x\ = —1. In projectable HL gravity, the gravitational degrees of 
freedom are the lapse function N(t), the shift function Ni(x,t ) and the 3-dimensional space 
metric gij(x,t), which appear in the ADM form of the space-time metric 

ds 2 = — c 2 N 2 dt 2 + gij [dx l + N l d£) (dad + N^dt) , (1) 

where c is the speed of light, with dimension [c] = 2. The mass dimensions of metric 
components are [N] = 0, [A,] = 2, and [g t j\ = 0. 

The covariant version of HL gravity involves the auxiliary gauge fields A(x,t) and v(x,t) 
whose role is to impose a constraint which eliminates the scalar degree of freedom of HL 
gravity, as explained below. The action is then given b}Q 

S = ^ J dtd 3 Xy/g^N - XK 2 - V + i/O ij (2K {j + V ; V,-//)| - AR j, (2) 

where 

Kij = 7^7 idij ~ ^i N j ~ v j N i } , hj = !, 2, 3 (3) 

Qij _ jv.i _ \ Rg ij 
2 y 

V = -c 2 R- ai R 2 - a 2 R ij Rij - ^iR 3 -/3 2 RR ij R lj 
-P 3 R? i R*Ri e - p 4 RV 2 R - p 5 V i R jk V i R? k . 

and the different tensors correspond to the 3-dimensional metric g t j . Note that [A] = 4, 
[v\ = 1, [e] =0, and the potential V includes all the renormalisable operators even under 
parity. The dimensions of the various terms in the Lagrangian are 

[fl] = 2, [B 2 ] = 4. [B s ] = [Afi 2 ] = 6, [c 2 ] = 4, [a,] = 2, [ft] = 0 . (4) 

Note that the term R l ^ kl Rij k i does not appear, as the Weyl tensor in three dimensions auto¬ 
matically vanishes. 

The action (J2]) is invariant under the following transformations: 

• 3-dimensional diffeomorphism 

5t = fit) (5) 

8x* = ?{t,x) 

&9ij ^j£>i T f, &k9ij T f 9ij 

SN t = d i £ k N k + Z k d k N i + &g ij + fN i + fN i 

5N = f k d k N + fN + fN 

SA = t k d k A + fA + fA 

1 It has been shown in pf)] that the present extension to the original HL gravity is valid for any A. 
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U( 1) symmetry 


S a N = 0 

&aQij 0 

5 a Ni = NVia 

S a A = a — N l Via 

dais = a ( 6 ) 

where a is an arbitrary spacetime function. 

In what follows, we shall need to distinguish tensors contracted with g tJ from those contracted 
with the flat-space metric S,j . For two vectors ly, w*, we will write the former as v l Wi and the 
latter as which leads to i+wy = —v 0 w 0 + v‘ l w i , where 

v l Wi = ViWjg lJ = (g~ 1 y :, ViWj = (l - ehij + e 2 h ik h jk ) ViWj + C>(e 3 ) . (7) 

Also, we denote v 2 = igi’i, we use <9* for the flat 3-space derivative, and d 2 for the flat 3-space 
Laplacian. 

2.2 Gauge fixing and £/(l)-symmetry constraints 

An interesting way of counting the dof of the theory is given in [18], where the Authors 
use a Hamiltonian description of gauge theories. They show that the number of primary 
constraints to take into account is the number of gauge functions plus the number of their 
time derivatives, in the situation where these gauge functions depend on space and time. 
This is because gauge functions and their time derivatives must be considered independent, 
when defining a boundary condition for the evolution of gauge fields. In our case, we have 10 
independent metric components (N, Ni, gij) and we see from the gauge transformations (l5il 
that the functions £* count twice since they appear with their time derivative, while / counts 
once only because it depends on t only. The total number of dof is therefore 10—(2 x 3+1) = 3. 
We are then left with 3 dof, consistently with the appearance of the Goldstone boson after 
breaking of the 4-dimensional diffeomorphism to 3-dimensional diffeomorphism. 

The metric fluctuations are defined by g l3 = Sij + ehij and, choosing the synchronous 
gauge, where N — 1 and = 0, we decompose h i3 as 

1 d d 1 

hij = Hij + diVj + djVi + (—Sij -+ -Sijh, ( 8 ) 

where is transverse traceless, Vi is transverse and h =tr {hij}. 

One can easily see that the variation of the action ((2]) with respect to A leads to R = 0, 
which is a condition we will impose in the path integral over graviton dof. This condition 
should be satisfied at the linear order in metric fluctuations, since we consider only quadratic 
terms in the gravity action. We obtain then 

0 = R = -e-d 2 (B + h) + 0(e 2 ) , (9) 

3 

which, together with boundary conditions h( oo) = B( oo) = 0 leads to h — —B everywhere. 
We can re-write our expansion of h,j as 

r\ r\ 

hij — Hij + diVj + djVi -|— ~^~h- (10) 
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Finally, We fix the U( 1) symmetry by setting u — 0, and we also note that ghosts 
decouple from matter at one loop, since the corresponding action is cubic in fluctuations of 
ghosts/gravitons. 

3 Effective dispersion relation for matter fields 

3.1 Coupling to matter 

We now wish to couple the gravity sector to classical matter, including a complex scalar 
held 0 and an Abelian gauge held Ap. In [19j a generic action is derived for the coupling 
of Covariant HL gravity to a scalar held, but we restrict ourselves to the case which recovers 
Lorentz symmetry in the IR 

Sscalar = ~ J d 3 xdt ^(-00* + (? (f' 3 8^8^) , (11) 

where c is the speed of light, with mass dimension 2. 

The coupling to an Abelian gauge held A M is described by the action 

Sphoton = ~ J d 3 xdt y/g(-2g 13 FoiF 0j + c 2 g lk g 31 F i:j F kl ) , (12) 

and we need not worry about gauge hxing for the Abelian held Ap, since it is considered an 
external source. We wish to calculate effective dispersion relations for the matter helds 0 
and Ap. As these helds couple to gravity only through their hrst derivatives, we may treat 
those derivatives as constant external helds: 

• 0 = 0o exp (ik^Xp) leads to 8^8^* = (fc) 2 0 

• Ai = A° sin (k^Xp) leads to F 2 - = 2 (k) 2 (A 0 ) 2 — 2 (k ■ A 0 ) 2 + 0(k 4 ). 

The (anisotropic) effective action will be given by an expression of the form 

S e JJ dar = j d 3 xdt(( 1 + a)00* - (1 + b)c 2 8i(j)di(j)*) , (13) 

with a corresponding dispersion relation 

(1 + a)kl = (1 + b)c 2 k 2 , (14) 

where we assuming a,b 1. If we note v 2 the product (phase velocity x group velocity), 
we have then 

v 2 

-l-l = b-a + 0(e 3 ) (15) 

c z 

Similarly, for the photon we will obtain an effective action of the form 

S;lL, = -\j d 3 xdt(- 2(1 + a')F K F m + c 2 (l + bf)F ij F ij ) , (16) 

with the corresponding correction 

^-1 = b'- a ' + Q{e 3 ) . (17) 

c z 

As noted in [13], the measurable quantity which violates Lorentz symmetry is the difference 

5v 2 = \v\ — v^\ = c 2 \b' — b — a' + a\ . (18) 
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3.2 Features of the one-loop integration 

Taking into account the above gauge fixing conditions and the constraint R = 0 in the 
path integral, the gravity sector expanded to second order in the metric fluctuations reads 

s gravity = ~ j d 3 xdt H tj (x) ( d 2 + 4 a 2 d 4 - 4/3 5 <9 6 ) H l3 

+ 2 Vjd 2 djVj + (1 - A )hd 2 h + 0(e) (19) 

Note that only the terms from the potential with no powers of R survive, and only those 
with two powers of R l3 contribute at this level, as R, l3 is 0(e). We can now see that we have 
only two propagating degrees of freedom (the two polarizations of H, 3 ), as expected. We may 
treat the cases of photons and scalars separately at the order at which we are working, as 
they have no interactions other than through gravity. 


For the scalar we have 


S, 


scalar 


= j d 3 xdt <( o + < -c 2 ( — -HijHii — -diVjdiVj — -h 2 ) d t 
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+ 0(e 3 ) \ (20) 


. ^ 

where is proportional to —\(j>\ + c 2 dtcpdicj)*, the usual relativistic Lagrangian density 
for the scalar. Following [13j, because of isotropy in space coordinates, terms of the form 
T l idi<fidj(j)* that are quadratic in the graviton fields have been replaced by (1/3)T/<9j</>ch0*. 
Terms that mix different components of the graviton cannot contribute to our corrections at 
0(e 2 ) and so are neglected here. Terms linear in the graviton fields do not contribute here 
either, as we treat the scalars as classical external fields; for the same reason, we are able to 
integrate the quadratic terms by parts. 

Similarly, for the photon we have 

Sphoton = \j d 3 xdt |£ ao + eV ( - l -R l3 H %3 - ^diVjdiVj - h 2 ^j F kl F kl + 0(e 3 ) j , (21) 

where isotropy has been used and Fao is proportional to the relativistic Lagrangian for 
photons. 


3.3 Integration 

As the components of the graviton do not mix, we may consider their contributions 
separately. 

3.3.1 Spin 2 


As we seek the difference between the modifications to the space and time components, 
we can neglect in the action the term q. We therefore need only to consider the action 


5 _ 1 , 

^ scalar ^ / 

d 3 xdt Hi 3 (x) | 

d 2 + 4a 2 d 4 

- 4(3 5 d 6 + ^d k (bd k ^ 

1 . 

(22) 

for the scalar and 

5 - 1 

^photon ^ i 

[ d 3 xdt Hij(x ) 

(d' 2 + 4 a 2 d 4 

-4 fed*+ 2 -^c 2 - A F kl F k ) 

1 H.iW ■ 

(23) 
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for the photon. 


Scalar field 

We denote by D(q) the Fourier transform of d 2 + 4ck 2 <9 4 — 4/3 5 <9 6 , and we have 


S H = - t 


d 3 pdpo d 3 qdq 0 
4,/ (2vr) 4 (2vr) 4 


%(p) [D(q) + %-c 2 d k 


6(p + q) Hij(q) . (24) 


Integrating over the two components of //, this gives a contribution to the partition function 
of 

(25) 


4e 2 


Det <J (D(q) + c 2 <9 fc 0<9 fc 0*)5(p + g) 


In ( D(q)S(p + q)) + -^-c 2 <9 fc 


= exp Tr 
= exp j-^c 2 <9 fc 0<9 fc </>*5(O) 


d 3 pdp 0 


* £ 1 (<?)<5(p + g) 

1 


(27t) 4 + 4« 2 p 4 + Af3 5 p e 


+ 


where 5(0) is a constant global volume factor, and dots represent held-independent terms. 
The integral in the above leads to: 


i = -m^-c 2 d k 


d 3 pdpo 


(27r) 4 -pi + 4a 2 p 4 + 4 fi 5 p e 


Af2 - OHT 

= - 5 Wl(2^ C x — 


d 3 p 


p- \/ r 'j - 


(26) 


and is logarithmically divergent. We regularize this integral by dimensional regularization, 
as it respects the symmetries of the theory 


L = -5(0) 


4e 2 


3(2tt)‘ 


-Jdkfidkfi* x 


—in f (3 — a)n 2 p 2 a dp 

2 J T(^)p 2 y/a 2 + P 5 p 


-P 


= *5(0)— jc 2 d k (j)d k (t) 

6n z V p 5 a 


* 1 n a r . 

+ hmte , 


(27) 


where p is an arbitrary mass scale. 

After dividing by *5(0), in order to take into account the Wick rotation and the space time 
volume, the identification with the speed (TT5|) gives 


v 4 -1 = 


—€~ 




67T 2 V^5 a 


(28) 


Photon 

Comparing the coefficients of the relevant terms in the actions, one can see the effective 
change in velocity for the photon will be 1/2 times that of the scalar. This leads to 



Yln 2 \ffil a 


(29) 


Note that both these results are sub-luminal, as might intuitively be expected: small fluctu¬ 
ations of the spatial metric about flat space should generally lead to a “longer” path between 
two points. 
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3.3.2 Spin 1 


The spin 1 terms are similar to the spin 2 terms shown above, with 2in the place 
of HijHij. However, the calculation leads to an integral of the form 


/ d 3 pdp 0 1 

WST 


(30) 


Unlike the previous case, the p integral here is the integral of a polynomial, which can be 
formally taken to vanish under dimensional regularization [[2D]. The vanishing or finiteness of 
a regularised integral which otherwise would naively be divergent is explained pedagogically 
in [2T|: in the regularised integral, divergences associated to different regions of the domain 
of integration cancel each other, such that the integral is finite when the regulator is removed. 


3.3.3 Spin 0 

The coefficients of the relevant terms and thus the modifications to the velocities coming 
from the spin 0 component of the graviton, h, are equal for the scalar and photonH Hence 
the final result below does not depend on the parameter A, which only appears in the spin 0 
kinetic term. 


4 Analysis and conclusions 

From the results (1281) and fl29l) . the total measurable difference in the squared velocities 


is 


— ^ 2^,2 ..a 

Sv 2 = v\ — v\ = „ -r= -- V finite . 


A I2n 2 y/Jk a 
We define the beta function, in the limit a —> 0, by 

d(Sv 2 


P = /*- 


e 2 c 2 


dp I2n 2 y/fk ’ 
and we can then write, for some mass scale p 0 , 

O 

p 


(31) 


(32) 


-e 2 c 2 

2 - 6 In 


Sv 2 = 


12vr 2 v^ ho 


In order to chose /i 0 , we repeat the calculation of the integral 
momentum cut-off A to obtain 


(33) 

regularised by the high 


67T 2 



(34) 


2 This is unsurprising, see the similar calculation for regular Horava Lifshitz gravity in Ha. the field the 
authors call a vanishes here. 













such that the identification with the speed (fT5l) gives 


1 


-e 2 A 2 
127T 2 V^5 «2 


( 35 ) 


and suggests po = \f &2 as a natural choice. 

The result (135]) is obtained in anisotropic Minkowski space time, and we rescale the time 
coordinate as t —» t Mjj L , where Mhl is a scale characteristic of HL gravity, below which the 
classical model can be considered relativistic. Speeds are then rescaled as v —> v Mjj 2 L , and 
we set the speed of light in isotropic Minkowski space time to 1 (we therefore identify the 
dimensionful quantity c = Mjj L ). The coupling constant appearing in the action (J2J) is then 
e = Mjjl/Mpi , where Mpi is the Planck mass, and the measurable deviation from Lorentz 
symmetry in the effective theory is 


Sv 2 = \vl-v\\ 


M 2 H l 


In — 

24vr 2 M 2 zV ^ 11 «2 


(36) 


which should be less than about ICG 20 [22]. 

In order to get an idea of the order of magnitude for Mhl, one can set the different 
parameters to natural values in the present context, which are (3$ ~ 1, a 2 — Mhl an d 
H — Mpi. This shows that one should respect M HL < 10 10 GeV for the result (1361) to satisfy 
the upper bounds on Lorentz-symmetry violation. This value is also obtained in [23], where 
non-relativistic corrections to matter kinetic terms are calculated in the framework of another 
4-dimensional diffeomorphism breaking gravity model, and where 10 10 GeV corresponds to 
the cut off above which the model is no longer valid. We note that this scale also corresponds 
to the Higgs potential instability [24], which could be avoided by taking into account curvature 
effects in the calculation of the Higgs potential [25], and it would be interesting to look for 
a stabilizing mechanism in the framework of non-relativistic gravity models. 

We comment here on the relevance of Lorentz-symmetry violation in the study of cos¬ 
mic rays, with energies necessarily lower than the Greisen-Zatsepin-Kuzmin cutoff Eqzk — 
lQi9.6i±o.o3gV |i26], since the latter is of the order of the bound we find for Mhl■ Above this en¬ 
ergy, protons interact with the Cosmic Microwave Backgnxmd, producing pions which decay 
and generate showers observable on Earth. As noted in [27], resulting photon-induced show¬ 
ers would be highly sensitive to Lorentz-symmetry violation, and the observation of 10 19 eV 
photons would put strong bounds on the different Lorentz-symmetry violating parameters of 
the Standard Model Extension - SME -[28]. The observation of such photons would also help 
put bounds on parameters in the expression f[36j) . which could then be related to the SME. 
The bound Sv 2 < 10"~ 25 found in [27] from potential ultra-high-energy photons, is actually 
much smaller than the one considered above, and concerns electron/positrons created by 
these high energy photons in the presence of the Earth magnetic field. Taking into account 
this bound, and assuming /3 5 of order 1, we find the typical upper bound M H l < 10 6 GeV. 
Although this is still well above the current accessible energies at CERN, the observation 
of ~ 10 19 eV photons would definitely put much stronger constraints on Lorentz-symmetry 
violation. 

We also remark that supersymmetric models have been studied in the context of Lifshitz- 
type theories [29], where an interesting feature is that particles in the same supermultiplet 
see the same limiting speed. Furthermore, in the case of supersymmetric gauge theory, the 
limiting speed is the same for matter and gauge supermultiplets. From the phenomenological 
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point of view, it is suggested that the Lorentz-symmetry violation scale Mhl should be 
above the SUSY breaking scale, in order to avoid fine-tuning problems imposed by bounds 
on Lorentz-symmetry violation. This means that Mhl should be at least of the order of 10 
TeV, which is well below the bound we find here. 

We make a final comment regarding the calculation carried out in |T3], for the original 
Horava-Lifshitz gravity and where a quadratic divergence is found. This stems not from the 
extra degree of freedom, but from treating the matter as quantum fields and then considering 
terms quartic in the matter fields, obtained from completing the square for the coupling 
terms linear in the graviton fields, that we neglected above. In terms of Feynman diagrams, 
treating matter as quantum fields consists in considering self energy graphs with an internal 
matter propagator. Our analysis shows that, as long as matter is classical, only logarithmic 
divergences arise. Studies similar to the present one are planned, involving higher-order 
space derivatives of matter fields, in order to generate dynamically Lifshitz-type theories in 
flat space time, as the ones analysed in HD El for example. 
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